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Abstract 

X) . 

rS^ ' We analyse the sum rules describing the action and energy in the colour fields around 

I glueballs, torelons and static potentials. 

(N , 

^ ■ 1 Introduction 

O 

■<!:;;j- ' In lattice gauge theory, it is possible to derive sum rules for the energy and action in the 

■ colour fields around states. The technique used in refQ — hereafter referred to as I — involves 
evaluating derivatives with respect to a parameter (/? for example) of the formal expression for 
a correlation involving an observable of intetest. For the Wilson gauge action, this yields exact 

■ relations between the /3-dependence of observables and the sum over a time-slice of the plaquette 
^ , expectation in the presence of that observable (for example see eqs. |l| and 23). These identities 



0^ 



Q I can be used as checks of numerical results — see ref||2[ for an application to the gluelump 

' state. They relate the variation with /3 to sums at one fixed value of (3. They can also be 

used to investigate the /3-dependence of lattice quantities: so leading to evaluation of the lattice 
^ ; beta-function. 

I A more powerful set of relations can be derived if the /3-derivatives can be re-expressed using 

renormalisation group invariance in terms of well known quantities. For glueballs and potentials, 
these were also presented in I. The main conclusion is that the combination of squared colour 
field strengths corresponding to the action (electric plus magnetic) is much larger than the 
combination corresponding to the energy (electric minus magnetic). This implies that the electric 
and magnetic field strengths are comparable. This conclusion has been a useful benchmark for 
models of non-perturbative QCD. Although this general conclusion was correct, the explicit 
results in I were in error and the correct expressions are given here. 

The lattice analysis of the field strengths depends on the scale at which these fields are 
probed. The results can be calculated reliably by perturbation theory for very short distance 
scales. However, for scales appropriate for non-perturbative states, perturbation theory on the 
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lattice in terms of the bare coupling is now known to be poorly convergent and it is worthwhile 
to re-assess the assumptions leading to these relations. 

Recently there has been a reawakening of interest in this area - partly from new accurate 
lattice results Q and partly because of the realisation [Q, ^ that the application of the sum 
rules to static potentials was wrong. 

Here we summarise the derivation of I and confirm the correction needed for the application 
to potentials. We extend the discussion to apply the sum rules to torelons and to analyse the 
transverse and longitudinal colour fields separately. This enables us to explore in detail the 
problem of the field energy in the potential between static sources. 



2 Glueballs 



As an example of the techniques to be used, we consider first a glueball state. Define M{(3) as 
the lattice observable glueball mass (in lattice units) which will depend on /3 the bare lattice 
coupling parameter, with /? = 2N /g^ ^ for the gauge sector of the SU(N) theory. 

Then, for the Wilson action, the identity was derived in I that 

^=<iiE°|i>-<oiE°io>=E°i-o (1) 

where □ is the plaquette action ■^Tr(l — Uo) which is summed over all (6L'^) plaquettes in one 
time slice. The subscript 1 — refers to the difference of this plaquette sum in a one glueball 
state (1) and in the vacuum (0). 

This identity can be used as it stands to check this observed plaquette difference with the left 
hand side obtained as a finite difference from lattice calculations of M{[3) at two nearby values 
of 13. A more powerful application comes from using the renormalisation group invariance to 
relate the /3-dependence of M to the /3- dependence of the lattice spacing a. Since M {(3(a)) /a 
is the physical continuum mass m as a — > 0, it must be independent of a. Hence 

^_dM{(3{a))/a _ M I dM d(3 
da a^ a df3 da 

Thus 

^ = E °i-o (3) 

ama 

Note that dp/dlna = — llA'"^/(127r^) to lowest order in perturbation theory for an SU(N) gauge 
theory. Thus the plaquette action is lowered in the glueball surroundings compared to the 
vacuum. 
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This is one of the prototype lattice action sum rules. It relates the plaquette action around a 
glueball to the mass of the glueball. It is exact provided that a non-perturbative determination 
of the lattice beta function is used. 

Further relations can be derived by splitting the lattice Wilson action into several terms 
with different coefficients. This analysis of asymmetric lattices was incorrect in I. In order to 
establish clearly the correct expressions, here we use a more direct method of derivation which 
also has the advantage of being more general. 

Consider the general case where there are different coefficients for all 6 orientations of pla- 
quette: 

^ E - E (4) 

There will be four lattice spacings Cj in general. We shall need to evaluate the derivatives 
df5ij/dak. At the symmetry point where ai = a for all i = 1,4, these derivatives fall into two 
classes 



S \i k = i OT j and — - — — = U \i k ^ i or j (5) 



d In d In 

The generalisation of the identities derived in I are also needed: 

dM 



E(°*j)i-o (6) 



where the sum is again over one time slice. Then the renormalisation group invariance of the 
result obtained on such a lattice implies that 

d M{(3jk{ao,ai,a2, a3),-) _q 



dai ao 

where oq enters because the glueball correlation is conventionally determined in the time direc- 
tion. Because only this time direction is privileged in this case, at the symmmetry point, we 
have that 

□oj = and Dj-fc = for j,k^O (8) 
where the subscript 1-0 is implied hereon. 

Applying the renormalisation group invariance conditions of eq. for i = and for i ^ 
gives 

M = J2^Sat + wa, (9) 

= Y,i'^U + S)at + {U + 2S)as (10) 
Then combining eq. |^ and eq. |l^ yields 

M = ^2(5 + [/)(3nt + 3n,) (11) 
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which is the the same as eq. H provided we have the consistency condition 



2{S + U) = -^ (12) 
a ma 



Subtracting eqsM and then gives 



M = J2^{S -U){3at-3a,) (13) 

This latter equation is appropriate to the energy in the colour field around a glueball. In 
order to make it more useful, we need to estimate the combination of derivatives S — U. 

Consider the special case, as used by Karsch P, where a* = oq; oi = 02 = 03 = ag and 
Poi = Pu Pij = Ps where i,j > 0. The derivatives in this case can be related to S and U, at the 
symmetry point: 

^ = 5and^ = S + 2C/ (14) 
a m at o In Qg 

f/and-?^ = 25 + [/ (15) 



din at d In a^ 

The dependence of Ps and Pt on at and a^ coming from the weak coupling limit of the 
theory is that, where ^ = Qg/at, 

Pt = CiPias) + 2Ncs{0 + ••) and Ps = r\P{as) + 2Nct{0 + ••) (16) 
where at ^ = 1: = q = and Karsch obtains = 0.114 for N = 2 and = 0.2016 for = 3. 



Then substituting eqs. 16 into eqs. 14 and 15 gives the constraint 



4iV(4 + c^) = -^ (17) 
a in a 

This is the same constraint as found by Karsch from similar consistency arguments. 
Using these expressions gives 

S = -0 + 2Nc!s + \-^ , U = P-2Nc', (18) 
2 am a 



and thus 



S-U = -2p + ANc'^ + l-^ (19) 

2 am a 



This implies that, as /3 — >■ 00, the energy sum-rule (eq. O) becomes 



^ = E ^Pi^°s - 3nt) (20) 
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The expression of eq. ^ in I had a factor of 1 instead of 4/3, coming from an error in the 
evaluation of the weak couphng result for the dependence of the asymmetric /3's on the a's. Note 
that the naive continuum expression for the energy in the colour field would be obtained with a 
factor of 1. 

In principle S — U can be determined non-perturbatively by simulating a lattice with non- 
equal /3's and determining the ratio of the lattice spacings in the 4 directions from the glueball 
correlations in those directions. Accurate data do not exist at present, although an indirect 
method has been used in SU(2) and substantial corrections are found Q to the weak coupling 
results. This is not surprising since lattice perturbation theory in the bare coupling is now known 
to be poorly convergent. This non-perturbative evaluation Q gives values of {U — S)/{2(3) of 
0.66 at /? = 2.4 and 0.77 at (3 = 2.8 compared to the weak coupling values of 0.85 and 0.87 
respectively. It is amusing that these non-perturbative estimates are close to 0.75 which would 
give the naive energy relation: 

M = ^/3(3n,-3nt) (21) 

The gluonic vacuum in QCD is known to be polarisable. It behaves like a medium and can 
be assigned an effective dielectric constant. Thus it is not really surprising that the naive sum of 
the energy in the colour fields (ie T,(3{30g — 30^) ) does not agree exactly with the mass. Indeed 
the result will depend in QCD on the scale at which the field energy is evaluated. A sensible 
scale would be commensurate with the glueball mass - where a non-perturbative determination 
of S — U is needed and rough agreement is obtained between the apparent field energy and the 
mass. The weak coupling calculation (which shows that only 3/4 of the mass lies as apparent 
energy in the colour fields) implies a very short distance scale of energy determination - which 
will probe the vacuum polarisation in a different manner. 

It is worth emphasizing the basic result, which was obtained in I already, that the electric 
and magnetic field strengths are comparable. In detail, the departure from equality is correctly 
given by 

S at -U-2S 3 dp 33N , , 

- = — = ^\ — — f« 1 (22) 

B 2U + S 2f5d\na 12-k 47r ^ ' 

where the approximation used in estimating S and U is valid at large (3. 



3 Potentials and torelons 



Having calibrated the approach on the glueball, we consider string states. The potential between 
static quarks is the case of greatest practical interest. Another related situation is with a closed 
loop of colour flux encircling the periodic boundary condtions: the torelon. As discussed in I, 
there is some subtlety in principle in dealing with the self energy of the static quarks. For clarity 
of presentation, the torelon case is considered first since the derivation is more compact. 
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The torelon is a closed string of colour flux in the fundamental representation that encircles 
the periodic boundary conditions in the x-direction where there are R lattice spacings in this 
direction. Its energy is measured on a lattice by analysing correlations of closed Polyakov line 
operators at t = and t = T. The study of the large T-behaviour then gives the lattice 
observable E{R,(3). 



The analysis of I gives, where R is kept constant, 

dE{R) 



dp 



E°i-o (23) 



where '1' now refers to the plaquette expectation value between torelon states and the sum is 
again over all plaquettes in a time-slice. The renormalisation group analysis now needs to take 
account of the fact that r = Ra must be kept constant in taking the limit a ^ 0. So 



dE{R,(3{a))/a 



da 

Thus 



E RdE ld(3 dE 
a^ dR a da 8(3 



(24) 

R 



^<«) + al^-5^S°.-o (25) 



As pointed out by Dosch et al. this expression differs from that in I where the term with 
the derivative with respect to R was omitted. The net effect of that term is, for a confining 
potential, to increase the effective left hand side of the action sum rule by a factor of 2. 

We now apply the general consideration of 6 couplings Pij as above. The new feature is 
that the torelon correlator is extended in the x and t directions. Thus we need to distinguish 
the X (longitudinal: L) and y, z (transverse: P) spatial directions. The 4 independent types of 
plaquette have orientations tL, tP, LP, and PP: we label them as £l, £p, Bp, Bl respectively 
in a natural notation. Note that £ here is related to the difference of the plaquette value in 
the torelon state and in the vacuum and so is the difference of gauge invariant combinations of 
electric colour fields squared. Following the same steps as above we obtain three independent 
constraints from the invariance with respect to ao, ai and ap of -^E{RaL, (5ij{a};)) 

E = J2S£l + 2S£p + 2UBp + UBl (26) 

a TP 

R—, = yS£L + 2U£p + 2SBp + UBl (27) 
oR ^ 

Q = Y,U£l + {S + U)£p + {S + U)Bp + SBl (28) 
where the sum is over one time slice. 



Combining these equations we obtain 



a JP 

E + R— = Y,2{S + U){£l + 2£p + 2Bp + Bl) (29) 
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E + R^ = Y.^{S-U){£l-Bl) (30) 

E-R^ = Y.'^{S-U){8p-Bp) (31) 

It is also convenient to write down the combination corresponding naively to the total energy 
in the fields: 

E - 3^y^ =T.^(S- U){£l + 2£p - 2Bp - Bl) (32) 

Again the action sum-rule (eq. ^9|) agrees with the result obtained from a symmetric lattice 
(eq. |25| ) with the same relationship oi S + U to the beta function (eq. [T^ ) as for the glueball case. 
Thus, apart from the term with a derivative with respect to i?, the results for the total action 



(eq. 2£) and total energy (eq. p^) are similar in normalisation to the glueball case introduced 



above. 

Consider, for orientation, the case where the torelon energy is a sum of a string tension piece 
and a string fluctuation piece: 

E{R)=KR-f/R (33) 

where we expect / = 7r/3 (this behaviour of the torelon energy has been checked numerically 
recently [^). Then the sum rules become 

2KR = Y,2{S + U){£l + 2£p + 2Bp + Bl) (34) 

2KR = Y,2{S-U){£l-Bl) (35) 
-2f/R = Y,2{S-U){£p-Bp) (36) 

This shows that the transverse energy in the fields (where here we define energy as £ — B ) will 
be much smaller for large R than the longitudinal energy. Moreover it has the opposite sign. 
These sum rules provide an independent way to study the split of the total torelon energy into 
string tension and string fluctuation components. 

Consider the sum rule for the longitudinal energy with the weak coupling (/3 — > oo) value 
for S-U: 

\kR = Y,P{Bl-£l) (37) 

The right hand side is just the naive expression for the energy in the longitudinal fields. Thus 
we obtain one half of the expected semi-classical result of KR. This is somewhat surprising 
since the longitudinal colour flux is applied explicitly and in the semi-classical limit the energy 
should remain relatively unaffected by quantum corrections. However, the vacuum polarisation 
effects will be strong at the large energy scale (corresponding to /? ^ oo) used to evaluate the 
expression. 

The application to the potential between static sources follows the same steps as for the 
torelon. The difference is that R is now the spatial extent of a Wilson loop rather than the 
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spatial extent of the lattice itself. The main new feature is that there will also be a self-enegry 
contribution in the lattice observable energy E{R). This self-energy was discussed in I. As a ^ 
it becomes the dominant term in the energy but it is very localised spatially. Thus it is possible 
to separate it out — leaving just the same results as for the torelon discussion above. One way 
to remove the self energy contribution, in practice, is by taking the difference of expressions for 
two values of R when it cancels. 



The analogue of eq. |3^ for the potential energy E{R) between static sources at separation R 



is a sum of self-energy, Coulombic and string tension terms: 

E{R) = Vo-e/R + KR (38) 



After removing the self-energy part (Vq), the sum rules then become the same as eqs. 34, |35| and 
|3^ with / changed to e. Thus the same result applies that the transverse energy in the colour 
fields will be much smaller for large R than the longitudinal energy (where here we define energy 
as £ — B ). This is a new result. 



4 Conclusions 



We have studied the energy and action distribution in the colour fields around non-perturbative 
states. We use the semi-classical definition of these distributions and define the appropriate dif- 
ference of the plaquette combination evaluated in the non-perturbative state and in the vacuum. 
The lattice definition we use has energy and action determined by the plaquette. In principle it 
should be possible to define a quantity which characterises the energy and action in the colour 
fields of a state and which has a continuum limit. Such a defintion could be based, for example, 
on using a square Wilson loop of fixed physical size as a — > 0. This would probe the energy and 
action distributions at a fixed physical scale and so would give a result free of lattice artefacts. 
But, of course, such a definition would not satisfy the sum rules we have derived. 

The simplest result, which was obtained in I and which has been checked in numerical studies, 
is that the action in the colour fields is much larger than the energy. This follows because the 
derivatives of (3ij with respect to on an asymmetric lattice can be expressed in terms of two 
independent quantities S and U which can be estimated. The correct expression for the ratio, 
derived here, is 

Energy _ 3{S + U) ^ -3 d(i _ g'^ 33iV 
Action ~ {S-U) ^ Jp dlna ~ 4^ 24^ ^ 

for a glueball state (where the approximation used in estimating S and U is valid at large /?). 
For the interquark potential or for a torelon this ratio is approximately 3 times smaller still. 

The naive expectation is that the spatial sum of the energy density in the colour field around 
a state should equal the energy of the state itself. We correct the results given in I, and find 
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that, evaluated by the semi-classical expression, the sum of the energy density in the colour 
fields around a glueball (and torelon) is given by 3/4 (and 1/2 respectively) of the energy of the 
state times 2.0(3/(JJ — S). Thus no non-perturbative value for S — U can make both of these 
sums exactly equal to the energy. Evaluating the energy density sum at a large energy scale, 
we can use perturbation theory to obtain a field energy around a glueball (and torelon) which 
is 3/4 (and 1/2 respectively) of the energy of the state. These fractions are closer to one when 
a non-perturbative estimate at a lower energy scale is used to determine the field sums. The 
explanation for the departure of these relations from identity is most easily achieved by invoking 
the vacuum polarisation effects as producing an effective dielectric constant. Moreover this 
effective dielectric constant must be different in the glueball (spherical) and torelon (cylindrical) 
geometries. 

We also present sum rules for the longitudinal and transverse field energy in a string state: 
torelon or interquark potential. For the potential between static quarks, this implies that at 
large R the transverse energy in the colour fields will be much smaller than the longitudinal 
energy. It will be interesting to explore this in lattice studies. 
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